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ABSTRACT 

Centrifugal pressure of matter spiralling onto black holes, have long been 
known to produce standing or oscillating shocks. The post-shock disc puffs up 
in the form of a torus, which intercepts soft photons from the outer Keplerian 
disc and inverse Comptonizes them to produce hard photons. The post-shock 
region also produces jets. We study the interaction of both hard photons and 
soft photons, with rotating electron-positron jets. We show that hard photons 
from the post-shock torus are instrumental in acceleration of jets, while soft 
photons from the Keplerian disc is a better collimating agent. We also show 
that if the jets are launched closer to the black hole, relativistic and collimated 
jets are produced; if they are launched at larger distances both collimation 
and acceleration are less. We also show that if the shock location is relatively 
at larger distances from the black hole, collimation is better. 

Key words: Accretion, accretion discs - black hole physics - radiation mech- 
anism: general - radiative transfer - ISM: jets and outflows 



1 INTRODUCTION 



Jets in mi cro quasars as well as in qu asars shows rela tivistic terminal spe ed [e. g., GRS 



915 



105. 



Mirabel fc Rodriguez! Jl99J); 3C 273, 3C 345 



Zensus et al. 



( 19951 ): M87, 



Biretta 



(1993|)], though the actual acceleration process is an enigma. It is well accepted in the sci- 
entific community that jets around compact objects originate from the accretion disc ac- 
companying such compact objects. The study of interaction of radiation from the disc with 
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outflowing jets, is quite extensive. The radiation field produced by a disc, depends on the 
geometry as well as the dominant physical processes of the disc model. And hence the study 
of interaction of jets with radiation field produced by different disc models, will, in general, 
dra w di ff erent conclusions. 

Ickd (|l980f ) studied the effect of radiative accelera tion o f gas fl ow above a Keplerian 



disc, but did not consider the effect of radiation drag. iPiranl (|198^), while calculating the 
radiative acceleration of outflows about the rotation axis of thick accretion discs, found out 
that in order to accelerate outflows to 7 > 1.5 (where 7 is the bulk Lorentz factor), the 
funnels must b e short and steep, but such funnels are found to be unstable. In a very im- 



portant paper, 



Icke| (|1989l ) considered blobby jets about the a xis of sy mmetry of thin discs 



(Shakura & Sunyaev 1973; Novikov & Thorne 1973, hereafter INT73I ) and he obtained the 
'magic speed' of v m = 0.451c (c — the velocity of light), t>, m being the upper limit of terminal 



speed. (IFukuel ll996) extended this study for rotating flow above a thin disc and drew simi- 
lar conclusions, although for rotating flow, away from axis of symmetry the terminal speed 
was found out to be a little less than the magic speed of Icke. It was shown that rotating 
winds above thin disc, generally spreads as the radiation force along the azimuthal direction 
(generated because of disc rotation) increases the angular momentum of the flow, and are 
diffi cult to accelerate because of the presence of radiation drag term in the vertical direction 



(Taiima Rr. Fnkiif 



1996, 



1998) . Radiatively dri ven winds were also studied for radiation com- 



ing from a slim disc (jWatarai fc Fukuelll9 99). Inspite of different disc temperature profile 
from the thin disc, and also the inclusion of advection term, radiatively driven outflows from 
slim disc concluded, that these radiation fields will spread the outflows and will suppress 
the vertical motion. Later Fukue and his collaborators, achieved relativistic terminal speed 
and collimation for pair-plasma jets from a disc model which contains inner ADAF region 



(non- luminous) and outer slim disc (luminous), for non-rotat ing black hole (IF 



ikue et al 



20011 ). and repeated the same scheme for rotating black holes (jOrihara fc Fukueil2003r ). We 



are working in a different regime. We consider the TCAF ( Two C omponent Accretion Flow) 



CT95I ). which consists of less lu- 



disc model (see, Chakrabarti & Titarchuk 1995, hereafter 
minous out er Keplerian disc and more lum i nous post-shock torus in the hard spectral-state 



of the disc. 



Chattopadhvav fc Chakrabartil (j2000i ) calculated the radiative force on optically 



thin j ets in and around the axis of symmetr y for radiations coming from the post-shock 



torus. 



Chattopadhvav &: Chakrabartil (|2002a| ) investigated the issue of radiative accelera- 



tion of normal plasma jets, by radiations only from the post-shock torus of TCAF disc, 
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and concl uded that normal-plasma jets are indee d accelerated to mildly relativistic termi- 
nal speed. IChattopadhvav fc Chakrabartil (j2002b[ ) also reported that hard radiations from 
the post-shock tor us do not impose any upper limit for terminal speed. Chattopadhyay et 
al. (2004, hereafter P aper- J ), solved the equations of photo-hydrodynamics fo r on- axis pair 
dominated jets for radiations coming from the whole of TCAF disc model, 
eluded that terminal speed t?~0.9c is easily achieved if the shock-location in accretion x s 
few xl0r 9 (r g =Schwarzschild radius), a nd the lu minosity of the post-shock torus (=L C ) is 



con- 



few x 10% the Eddington luminosity. 



Paper-J also showed that the special geometry of 



the post-shock torus ensures that, its contribution to various radiative moments dominates 
over that due to Keplerian disc. Thus radiations from Keplerian disc (of luminosity Lk) has 
a very marginal effect on determining the terminal speed of jets. Having shown this, it is 
only natural to study how rotating jets behave in the presence of the radiation field of the 
disc. The azimuthal component of the radiative force generated by the disc rotation, will 
try to spin up the jet, while radiative drag force in the same direction will try to carry away 
the angular momentum from the jet. So it is particularly interesting to study which of these 
two phenomena wins. 



In this paper we study the interaction of radiations from the whole of TCAF disc and ro- 
tating jets. We show that the radiation from t he CENB OL (CENtrifugal pressure-dominated 



BOundary Layer; see Chakrabarti et al. 1996. 



CTKE96I ) is the chief accelerating agent, while 



radiation from Keplerian disc has greater degree of collimation property. Both collimation 
and acceleration are much higher when the jet materials are launched closer to the axis of 
symmetry. We show that, if radiative processes is the dominant accelerating agent then, we 
can get highly relativistic and collimated jets only in moderate hard states (i.e., L C ~L K ) 
and not in extreme hard states (L C ^>L K ). 



In the next Section, we present a brief account of the TCAF disc model. In §3, we present 
the model assumptions and the equations of motion. We give an account of the streamline 
coordinates which we use to solve the governing equations. We also compute ten independent 
moments of radiation field from such an disc. In §4, we present our solutions and finally in 
§5, we draw our conclusions. 
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Figure 1. Cross-sectional view of Two Component Accretion Disc Model. 



2 TCAF DISC 



A detailed acc ount of TCAF disc is given in Pape r- 1 [also see 



1997, hereafter 



C97 



Chattopadhvav et al. 



CT95 



CTKE96; Chakrabarti 



( 20031 )]. For the sake of completeness, let us now 
give a brief account of the TCAF disc. 

Inner boundary condition for matter accreting onto black holes are (i) supersonic and 
(ii) sub -Keplerian. Accretion topo logies for sub-Keplerian matter admits two X-type critical 



points (jLiang fc Thompson 1980), and supersonic m atter a 



point may undergo shock due to centrifugal pressure (jFukue 



ter crossing the outer critical 



1987 



Chakrabarti 



1989 



1990 



19961 ). The post-shock matter then enters the black hole through the inner critical point. 



Matter is slowed down in the immediate post-shock region and as entropy is generated, it 
makes the post-shock region hot. This hot, slowed down post-shock region puffs up in the 
form of a torus the CENBOL. 

Chakrabarti & Titarchuk (1995) for the first time proposed a disc model which contains 
both, Keplerian matter and sub-Keplerian matter. The Keplerian matter is of higher angular 
momentum and lower specific energy than the sub-Keplerian matter, and settles around the 
equatorial plane to form the Keplerian disc (hereafter KD). Sub-Keplerian matter flanks 
the cooler Keplerian disc from the top and bottom, sand wichi ng the Keplerian disc, and is 
known as the sub-Keplerian halo (hereafter SKH; also see. lC97l ). The SKH suffers a shock at 
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a few tens of Schwarzschild radii. The shock compresses the post-shock flow making it denser 
than the pre-shock sub-Keplerian flow. The hot post-shock flow evaporates the Keplerian 
disc and falls to the black hole as a single component, in other words the shock location (x s ) 
is the outer boundary of CENBOL and the inner boundary of the KD. A schematic diagram 
of such a disc structure is shown in Fig. (1), where KD is shown to be sandwiched by SKH, 
and the position of x s , as well as the central black hole are also shown. Although the SKH 
(pre-shock) is optically thin for the radiations from Keplerian disc, the post-shock torus is 



optica l 



CT95, 



y slim, in othe r words the optical depth of CENBOL is around unity (i.e., r~l; see 



CTKE20, for details). 



CT9^ showed that, if the Keplerian accretion rate (Mr) is higher than sub-Keplerian 



accretion rate (Mh), then it supplies more soft photons to cool down the CENBOL. This 
results in more power to the soft end of the accretion disc spectrum — a state known as 
the soft state. If on the other hand, Mh > Mr, then SKH supplies more hot electrons to 
the CENBOL than the soft photons supplied by the KD. The dearth of soft photons cannot 
cool down the CENBOL significantly. Thus CENBOL remaining puffed up and hot, can 
intercept a significant fraction of soft photons produced by KD, and inverse- Comptonize 
them to produce the hard power-law tail of the accretion disc spectrum — a state called 
hard state. This kind of hybrid disc structure is known as the Two Component Accretion 



Flow or the TCAF disc (1£ 



support (S mith et al. 



2001 



95; 



CTKEQf 



Ebisawa et al 



1996; 



C97l ). and has observational 



2002). 



Chakrabarti and his collaborators have also shown, the unbalanced thermal gradient 
force of CENBOL, in t he z direction, drives a part of the in falling matter along the axis o f 



symmetry to form jets ((Chakrabarti 



1998, 



1999; 



Das fc Chakrabarti 



1999; 



Das et al. 



2001). 



There are wide support that the jets a re indeed coming o ut from a region within 50 — 100 



Schwarzschild radius of the black hole (IJunor et al. 
are produced only in hard states [see 



Gallo et al 



1999). Similarly, it is believed that jets 



2003T ). and references therein]. Thus it 



is natural to study interaction of hard radiation from the CENBOL and the outflowing jets, 
with the particular interest of studying, whether momentum deposited to the jet material 
by these hard photons can accelerate them to ultra- relativist ic speeds. 

It is to be remembered that, we are not considering generation mechanism of jets self- 
consistently. Since hard radiations are expected to emerge out of CENBOL, the hard photons 
'look' directly into the jet vertically above and hence eventually deposit their momentum into 
the latter. Furthermore, radiation from a hot CENBOL is a likely source of pair production 
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and hence the possibility of radiative momentum deposition is likel y to be higher even for 
radiations from CENBOL hitting the outflow at an angle [see, e.g., (Y amasaki et al. 1 11999) 
for production mechanism of pairs from hot accretion flows]. Thus we consider radiative 
momentum deposition on pair dominated jets, which are generated above the inner surface 
of CENBOL, i.e., within the funnel like region. 



3 ASSUMPTIONS, GOVERNING EQUATIONS AND COMPUTATION OF 
RADIATIVE MOMENTS FROM TCAF DISC 

We ignore the curvature effects due to the presence of the central black hole mass. The 
metric is given by, 

ds 2 = c 2 dP - dr 2 - r 2 d(j) 2 - dz 2 , (1) 

where, r, <p, and z are the usual coordinates in cylindrical geometry and ds is the metric 
in four-space. The four-velocities are u^. The convention we follow is — the Greek indices 
signify all four components and the Latin indices represent only the spatial ones. The black 
hole is assumed to be non-rotating and hence the strong gr avity is taken care of by the 
so-called Paczyhski-Wiita potential ()Paczvhski fc Wiital ll98C0. 

In this paper, the generation mechanism of jets is not considered. We assume axis- 
symmetric, steady, rotating jets i.e., d/dt = d/d<f) = 0. The derivation of the equations of 
motion of radiation hydrodynamics for optically thin plasma, was invest igated by a number 



of workers. A detai 



ed ac count of such derivation has been presented by 



Mihalas fc Mihalas 



(J198J) and iFukud ([1996T 1 [references therein], and are not presented here. The equations of 
motion are; 



(e + p) u^— — - 4 



GM B r 
dx» ' R(R - r g ) 2 



— rivw 



dp 
dr 



dp ot_ 

dx 11 ^ mc 



(e + p) u> 



-u r u^ 



dp 
dxv 



+ P- 



mcr 



and 



. J u du z GM B z 



dz dx^ mc 



(2a) 
(26) 

(2c) 



In above equations, e, p and p are the internal energy, isotropic gas pressure and density 



measured in the co-moving frame of the fluid and R 



+ z 2 ) 1/2 . G, M B , a T , m and 



r g (2GMB/c 2 ) are the universal gravitational constant, the mass of the central black hole, 
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Thomson scattering cross-section, mass of the gas particles and Schwarzschild radius, respec- 
tively, signifies radiative contributions along respective components of the momentum 
balance equation. 

The radiative contributions are given by; 

pr 



m c 



Similarly, 



m c 



0~t 

m 



7 



— - 7 V£ - u ] P rj + u r ( 2^u j F j - u ] u k P jk 



if — l 2 u r e — u j p rj + u r (2^yu j f 



c 

'■' - u j u k p ik } 



0~t 

m 



7- 



p<t> 



- *pru+E - u j P^ + ru 4, (2-u j F j - u j u k P jk ^j 



,7 



1 f - rfrufe - v?p^ + ru 4 [ 2-v? f 3 - u j u k p jk 



and 



0~x 3 2 

m c 



0~t 

m 



7" 



- 7 V£ - u j P zj + u z ( 2-u j F j - u j u k P jk 



7 / 2 - 7 Ve - u j p zj + u z (2'ju j f j - vPvFp'^ 



(3a) 



(36) 



(3c) 



In Eqs. (3a-3c), E(r,(j),z), F l (r, </>, z), and P tj (r,(f),z) are the radiative energy density, 
three components of radiative flux and six components of radiative pressure tensors measured 
in observer frame, while e = —E, P = ^F l , and p lJ = sz-pv. Furthermore, 7(7 = u+) is 
the Lorentz factor. We assume the gas pressure to be negligible compared to the radiation 
pressure terms. 
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3. 0. 1 Streamlines 

In streamline coordinates (s, <ft, t), s-axis (t-axis) is parallel (perpendicular) to streamlines, 
where ds 2 = dr 2 + dz 2 , such that u l is zero. In Fig. (2) we draw a schem atic diagram 
of the streamline coordinates. The basic equations can then be written as (jFukud llf" 



Fukue e.t a,l. 



2001). 



Momentum Balance along the streamline: 



Bu s GM B rdr + zdz (r 2 vft) 2 dr a T ~ a , . 

u ^- = -^T5 ^ 1 + —^^ + — 3 > 4g 

as K(K — r g y ds r 6 ds cm 



Angular Momentum conservation: 



os cm 



And the streamline equation: 

GMb rdz — zdr A A dz , n „dz „^dr, 



+ ru^u^— + (f— - f z —) = 0, (4c) 



R(R — 2r g ) 2 ds ds ds ds' 

where, 



— — = if s - j 2 eu s - (p ss u s + rp s W) + u s {2 1 (f s u s + rf V) 
m c 

- (p ss u s u s + 2rp a *u'u 4> + r 2 p H u <t 'u 4 ')}. (Ad) 



In Eq. (4d), 



and, 



ds ds 



ss = rr fdA 2 + 2 p rz — — + ( —\ 2 
V ds I ds ds \ ds I 



ds ds 



Equations (4a-4c), can be simplified in terms of v and A and can be expressed in geometric 
units (r g = c = M B = 1), for simplicity we will keep the same symbols representing various 
quantities defined so far. 

We now define a three velocity measured by the static observer in geometrical units defined 
above, 
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while the angular-momentum and angular velocity are being defined as, 

A = -^ & tu = ^t. (5b) 

u- t u l 

We further define the 3-velocity measured by co-rotating observer, 

v 2 

v 2 = (5c) 

So that 7 2 = 1/(1 - v 2 - toX) = 1/{(1 - v 2 )(l - toX)} = 7 2 7 2 , u s = lv v, and = ("fX)/r 2 . 
Equations (4a-4c) are re-written as, 

dv JVi 



where, 



Ni = [ 7 J^ S - ^ v v£ - lv vV ss - 1 -V S * + 2 1 ( 1 2 u v 2 r s + llv v-F<t>) - 7 yp ss 
and 



L>! = 7%. (6c) 

£ = £ (<w) 

dz D 2 y ' 



N 2 = [7^-7 3 -^-7^^-7-^ + 7-{27(7^ s + 7-^)-(7y^ ss 

+ 2 7 7^-^ + 7 2 ^^)}]^-77yA^ + 77y^^ (6e) 



^ 2 = 77, (l + ^) « (6/) 



and, 

dr -r/[2R{R - l) 2 ] + 7 f + (7A) 2 /r 3 
~ -z/[2R(R- l) 2 ] +7/ 2 

Quantities defined in Eqs. (6b) and (6e) are: 



(6<?) 



r-r + r 

dz 
dz 



(7a) 
(7b) 
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Figure 3. Cartoon diagram of CENBOL. The outer KD and SKH are not shown. B(r, 0, z) is the field point and D(x,ip, 
y) is the source point on the inner surface of the CENBOL. The position of the black hole is at O. The inner surface of the 
CENBOL is assumed conical. The position of the black hole is O. ON = x s , and HH' = h s . LB'OD' = (j> & LD'OD" = d(j>. 
And dA is the differential area at D. Only the top half is shown. 



V s 



£ - e— 

dz 

rr dr dr rz dr z dz 

P -r-r + 2p 1~ + P T 

dz ds dz ds 



V s 



dz 



j>4 



P" 



,ds 
dz 



(7c) 
{Id) 
(7e) 
(7/) 



We now have to solve Eqs. (6a), (6d) and (6g), for given radiation field (e, /*, p u ), 
specified by disc parameters. It is to be noted that for motion along the axis A = and 



dr 
dz 



dX 

dz 



0, then Eq. (6a) reduces to Eq. (6) of 



Paper-] 



3.1 Computation of radiative moments from TCAF disc 

The radiation reaching each point within the funnel like region and the region above it, is 
coming from two parts of the disc, namely, the CENBOL and Keplerian disc, hence all the 
radiative moments should have both the contributions. 

In Fig. (3) a schematic diagram of the CENBOL is presented. The inner surface of the 
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CENBOL is assumed conical (described by, rotating OH about the axis of symmetry), and 
the outer surface of CENBOL is cylindrical (HH'). XOY describes the equatorial plane. 
O is the position of the black hole. B(r, 0, z) is the field point where the various radiative 
moments are to be calculated. D(x, <p, y) is the source point on the CENBOL inner surface. 
ON = OH' = x s is the location of the shock. And HH' = h s is the shock height. The local 
normal at D is m§r- The differential area about D is marked as dA. The shoc k height h s 



depends on the shock location, and is expressed as /i s ~0.6(x s — 1) [see. iPaper- 
By definition, the radiative moments (in natural units) at B are; 

E = -j ' Idti = - (J ' I c dKl c + J I K d£l K ) , (8a) 

Y = lJ Ilidn = \ (f Icicle + J J K V K dn^j , (86) 

and 

pv = I J m j dQ = - (J c i c r c i j c dn c + j i k IkI\c&Ik) • ( 8c ) 

In Eqs. (8a-8c), /, dfl, l l are the frequency integrated intensity from the disc, differential 

solid angle at B, and the direction cosines at B w.r.t D, respectively. Suffix C and K 

represent quantities linked to CENBOL and the KD respectively. The expressions of solid 

angles subtended at B from D is given by 

,^ (dA) cos LCDB x cosec^ dx dd> coslCDB , 

= BEfi = BW> ' (M) 

where, coslCDB = {M ^0j%0f 1 , and BD 2 = r 2 + x 2 + (z-y) 2 -2 r x cos0, CD = x secfl, 

BC 2 = r 2 + (z—x ta.n8—y) 2 and 6[= tan~ 1 (x s /h s )} is the semi-vertical angle of the CENBOL 

inner surface. Similarly it is easy to find from, Fig. (2), the direction cosines are given by, 

V c = ( r _ x cos0) /BD; 1% = -(x sin0) /BD; and l* c = (z - y)/BD. 

The frequency averaged CENBOL intensity is given by, 

j _ Icq _ L c 

C ~ (l + z md y ~ 7rA(l + z ved r 1 6> 

where, Lc and A is the total CENBOL luminosity and total CENBOL the 

red-shift factor taken up to the first order of the disc velocity. The inner edge of CENBOL 

i.e., also the TCAF disc is taken up to Xi n = 2r g , within which the general relati vistic eff ect 



Paper-] . 



has to be considered. Icq is assumed uniform for simplicity, as was explained in 

We are not solving the disc-equations simultaneously, so we, in principle, cannot consider 
the Doppler shift of the photons coming out of CENBOL. Nonetheless, not considering 
the disc motion, robs us of a vital element of physics. The rotational velocity of the disc 
generates the 4> component of the radiative flux. So we make an estimate of the post-shock 
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Figure 4. (a) Cartoon diagram of CENBOL. The outer Keplerian disc and sub-Keplerian halo are not shown. B(r, 0, z) is 
the field point and D(x,<f>, y) is the source point on the inner surface of the CENBOL. The inner surface of the CENBOL is 
assumed conical. The position of the black hole is also shown. ON = x s , and HH' = h s . LB' OD' = <j> & LD'OD" = d<j>. And 
dA is the differential area, (b) Cartoon diagram of CENBOL & Keplerian disc (KD) geometry. The SKH is not shown. B(r, 0, 
z) is the field point and Dfc(xfc,<j), 0) is the source point on the KD. The CENBOL is represented by the cylinder. The black 
spot represents O the position of the black hole. OH' = x s , and HH' = h s . And dAn is the differential area at Dk (shaded). 



disc motion. At the shock, the in-falling matter is virtually stopped. We solve the geodesic 
equation, starting from x s with a very small velocity (~0.01c). We assume the solution 
(u s t) [Appendix A] to be the 3-velocity of matter along the surface of the CENBOL. The 
angular momentum of th e post shock matter is sub-Keplerian and almost constant [e.g., 



Chakrabartil ()199 6)]: as the infall time-scale is much smaller than the viscous time- 



scale. So it is assumed to be constant, and is also the initial specific angular momentum 

(Aj n = 1.7 in the geometrical units defined above) of the jet. Under such assumptions, the 

radial, azimuthal and axial 3-velocity of the matter on the CENBOL surface is given by 

(expressed in dimensionless units described above), u x = u st sm9, u§ = (1 — l/x) 1 ^ 2 (Xi n /x) 

and u y = u st cos9. Therefore (1 + z re d) = (1 — Uil 1 ), where 

• (u x cos0 — tt0 sin0) (r — axos0) {u x sin0 + cos0) (x sin0) u y (z — y) 
Ud = BD BD + BD ■ {8f) 

In Fig. (4a-b), we represent a cartoon diagram of CENBOL and Keplerian disc (KD). In 

Fig. (4a), D k (xk,4>i 0) is the source point on the KD, and as in previous figure, B(r, 0,z) 

is field point where the various moment are computed. In Fig. (4b), only one half of the 

CENBOL/KD geometry above the equatorial plane is shown. B is still the field point, but 

D k (xk, 4>f, 0) is the limit, up to which B can see the annulus on KD defined by radius (%), 

in the positive <fi direction. So 0/ is the limit of integration for (f). From Fig. (4a), one can 
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easily find out, 

, n dA K /nnn zx K d<f)dx K 

^ = BD% C0SlCDB = BD\ ' ^ 
where, BD 2 K = r 2 + z 2 + x\ — 2r xk cos0, and the direction cosines are given by, 
V K — (r — a;x cos<p) / BD K ; l K = —(xk sin0) / BD K ; and = z/BD K . 
The frequency averaged KD intensity (see, NT73) is given by, 

Iko (3GM B M K )/(8n 2 r 3 g ) (x K 3 - ^x K 7 ' 2 ) 



4 



(8h) 



Now, the Doppler shift term is given by 1 + £ re d = 1 — W;X K , where, 



H _ . .r-x K cos(j) _ x^sin0 
CC7j/ x = -u K sm0 — u K COS0— — — . (8z) 

In the above equation -u^- = ^1/{2{x k — 1)} — the Keplerian velocity around a non- 
rotating black hole (in geometrical units). 
Shadow effect of CENBOL on the jets: 

As the jets are produced within the funnel like region of the TCAF disc, up to certain z 
the radiation from KD to the jet material, is blocked by the presence of CENBOL. It is easy 
to find from Figs. (3-4), that for material at a particular r (< x s , as jets are produced from 
the CENBOL region), the radiations from KD is completely blocked for z^h s (x — r)/(x — 
x s ). But even for z > h s (x — r)/(x — x s ), the jet cannot 'see' the entire Keplerian disc. 
From Fig. (4b), it is clear that the jet material at B cannot see, the whole of the annular 
area defined by radius xk, as part of it is blocked by the CENBOL. If for a particular r, 
h s (x — r)/(x — x s )^z^h s (xK + r)/(xK — x s ), then, from Fig. (4b), one can find a expression 
for /; 

where, B'D 2 K = r 2 + x\ — 2rx K cos<pf, and FD K = (h s B'D)/r. <pf is to be computed 
numerically from Eq. (8j). For, z > h s (x K + r)/(xx — x s ), 0/ = 2-k. Eqs. (8a-8c) are 
integrated, with the help of Eqs. (8d-8j), to get the expressions of various radiative moments 
of a TCAF disc. 

Let us now multiply, ar/m, with Eqs. (8a-8c), to get, 



£ — £ko 



z[x 



K 2 - Vsx^ w 



(r 2 + z 2 + x\ - 2rx K cos(p'fl 2 (l + £ red ) A 



[** z{x K 2 - y/hr^w , 

J 4 >f(r 2 + z 2 + x 2 K -2rx K cos ( j ) ')y 2 {l + t red )* 1 K 
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r 2 - (x 2 + y 2 f/ 2 cosiCDBd(t)'dx 



/x s rZTl 
x in J 



+ z 2 + x 2 - 2raxos0') 3/2 (i + z red ) 4 
= e K oE K (r, z, x s , x Q ) + e C oE c (r, z, x s , x ) 

= e K + e c (9a) 

p r*» r*f z[x- 2 - ^ 5/2 ) P K d(j>' 

7 " JK0 J Xs [ Jo (r 2 + z 2 + x 2 K -2rx K cos<f ) ') 3 / 2 (l + ^ ed ^ 

[i* z[x~ 2 - y/3^ 572 ) Z^' 

7 (r 2 + ^ 2 + x\ - 2r W os0') 3/2 (l + treaT % 
(x 2 + y^cosLCDBVdtfdx 

iC °l xj (r 2 + 2 2 + X 2 - 2^0080') 3/2 (1 + ^red) 4 

fmF l K {r; z, x s , x Q ) + fcoF l c (r, z, x s ,x Q ) 

fx + fb (96) 



+ 



Pko 



rx r<pj 

5 / U 

J x s J 



(r 2 + ^ 2 + ^ - 2r W os0') 3/2 (l + ^ed) 4 
z{x~ 2 - V3^ 5/2 ) & 

^ (r 2 + Z 2 + X 2 K _ 2rax COS0') 3/2 (l + Cred^ % 



xi n J 



2 - (x 2 + y 2 )^ 2 coslCDB l l Vdfidx 



' (r 2 + z 2 + a; 2 -2ra:cos0 / ) 3 / 2 (l + ^ e d) 4 
= p K0 P%(r, z, x s , x Q ) + pcoPc( r , z, x s , x ) 

= pl + Pc (9c) 

In Eqs. (9a-9c), the space-dependent part of e, / l , and p lJ are expressed as E 1 , and 
P u . Suffix 'X' and 'C signify Keplerian and CENBOL contributions. If the moments are 
expressed in dimensionless units then the constants in Eqs. (9a-9c), are given by, 

1.3x10 38 4^t , Q ,, 
eco = fco = Pco = 2ncmAGMG (M) 

and 

, 4.32xl0 17 m fc cr T c 

where, 4 = L c /L Edd (L c — CENBOL luminosity & L Edd — Eddington luminosity), 
rhk = Mk/Mecm (Mk — Keplerian accretion rate, ME dd — Eddington Accretion rate), and 
A is the CENBOL surface. 

For simplicity, we will not compute the shock location x s or the the CENBOL luminosity 
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(Lq) - instead, we w ill supply them as free parameters. They can be easily computed from ac- 



198 



9 



1990; 



CT95I 



DaS e-t al 



2001 



Chattoparihvav et al 



cretion parameters (jChakrabarti 

2nd). 

To obtained all the components of radiation field, we supply the following disc-parameters, 
(a) the inner radius of the CENBOL x- m {= 2r g , as explained in §3.1 ), (b) the shock loca- 
tion x s , (c) the CENBOL luminosity 4 (in units of L Edd ), (d) the Keplerian accretion rate 
rh k . 



The expression of Keplerian luminosity was given in 



J K 



27lI K0 27lXKdXK = -j-TTlfc 



Paper- II and is; 
3~ 



^Edd 



-fc-^Edd 



(10) 



It is to be remembered, that as jets are only observed in intermediate to hard spectral 
states of the accretion disc, so we will constrain our analysis in the domain 4/4~l to 

4/4 < i. 



3.2 The components of radiation field 

In Fig. (5), we show the contour plots of the space-dependent part of various moments of 
radiation due to the CENBOL. The shock location is x s = 20r g . The space variation of 
various moments due to the CENBOL is strongest within the funnel like region of the disc, 
so we only plot them within and above the funnel like region. Further as the jets are likely to 
be produced in this region, so this is the region that matters for our purpose. Contour plots 
of various moments are Ec(r,z) [Fig (5a); max. value: 50.55], F£,(r, z) [Fig. (5b); max/min 
value: 1.7/-5.42], Fq(t,z) [Fig. (5c); max. value: 28.94], F^(r,z) [Fig. (5d); max/min values: 
5.2/-1.19], P£ r (r,z) [Fig. (5e); max. value: 14.62], P r ^{r,z) [Fig. (5f); max/min values: 0.98/- 
1.66], pQ Z (r,z) [Fig. (5g); max/min values: 0.67/-1.39], P^(r,z) [Fig. (5h); max value: 29.8], 
pQ Z (r,z) [Fig. (5i); max value: 4.4] and P^ z (r,z) [Fig. (5j); max value: 3.88]. 

There are few features that are to be noted from the various moments of radiation, 
computed from the CENBOL: 

(i) The radiation field is highly anisotropic within the funnel like region of the CENBOL. 

(ii) All the moments maximizes and have very sharp gradients around the inner edge Xi n . 

(iii) Within the funnel {i.e., for r < x s & z ^ h s ) Fq < 0, but at moderate values of r & 
z away from the black hole, F r c > 0. 

(iv) F^ < very close to black hole. 

(v) F$ > for r > & z > 0. 





(j) 



Figure 6. The spacc-dcpcndcnt part of the ten independent components of radiation field due to the CENBOL. (a) Ex(r, z), 
(b) F- K (r,z), (c) F*(r,z), (d) F£.(r,z), (e) P r K r (r,z), (f) P r /(r,z), (g) P?(r,z), (h) P**(r,z), (i) P**(r,z), (j) /»*(r,*)- The 
disc parameters arc x s = 20r 9 & x = 500r 9 . 
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(vi) At r = 0, F$ = F c = 0. 

(vii) At r^x in k z^small, F^ > F c > F c . 

(viii) At r~Xj„ & z— >• large, Fq^hFq but less than F c . 

(ix) P^ is the most dominant of all the moments. 

(x) P r c r < & p r c r(p}f > ) > p* z . 

(xi) p^p- 

(xii) For x s ~10 — 20 r n fc R— »100rn, the radiation field due to CENBOL approaches that 



due to a point source (|Paper-II ) . 



It is evident that, as z becomes small, Fq > F c > Fq, still higher values of Eq & Pp 
will ensure that the A gained by Fq will be less than that reduced by the drag terms of 
Eq. (6d). It is also evident that as F c changes from < to > 0, as one goes away from the 
black hole and the axis of symmetry, which means closer to the axis and the black hole, 
radiation from CENBOL would push the jet material towards the axis and further away it 
will tend to spread the jet. This doesnot mean that, more luminous the CENBOL, more is 
the spreading, since radiation drag along r-direction will limit the radial expansion of jet. 

In Fig. (6), the contour maps of various moments of the radiation field due to the outer 
Keplerian disc is plotted. The shock location is same as Fig. (5). In contrast to the contri- 
bution due to the CENBOL, the Keplerian contribution is zero within the domain r^x s , 
z^h s (x Q — r)/(x a — x s ). The anisotropic nature of the radiation field extends to a region 
at much larger distances above the funnel like region of the CENBOL. The contour map of 
various moments from the Keplerian disc are, Ek{t,z) [Fig. (6a); max. value: 2.23xl0 -5 ], 
F T K (r,z) [Fig. (6b); max/min values: 1.77x l(r 6 /-1.12x 10~ 6 ], F&(r,z) [Fig. (6c); max. 
value: 8.62xl0~ 7 ], F^(r,z) [Fig. (6d); max. value: 1.7xl0~ 5 ], P^(r,z) [Fig. (6e); max. value: 
5.75xl(T 6 ], P$(r,z) [Fig. (6f); max/min values: 9.49xl(r 8 / - 1.37xl(T 8 ], P r K z (r,z) [Fig. 
(6g); max/min values: 1.25xKr 6 /-8.6xl(T 7 ], Pjffaz) [Fig. (6h); max. value: 5.75xl(T 6 ], 
Pp(r,z) [Fig. (6i); max value: 5.9xlCT 7 ] and P^ z (r,z) [Fig. (6j); max. value: 1.4xl0~ 5 ]. 

As in the previous figure, there are few features to be noted in Fig. (6) too. They are: 

(i) Each of the moments of the radiation field due to the KD is few orders of magnitude 
less than the corresponding ones due to CENBOL. 

(ii) Max(Pjf) > max(P^) > max(P^). That means spreading of jets will be less. 

(iii) The gradients of moments from the KD are lesser compared to that from the CEN- 
BOL. 
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(iv) F K > 0. 



(v) F T K < in a larger domain =>- pushing the jet materials towards axis in a larger 
domain above the funnel of the CENBOL, and the angular momentum gained will also be 
less as (Fx) is smallest of the three components of flux. 



(vi) P K r rzP K *fxP K z . 



It is quite evident that in the hard state, Keplerian contribution to the radiative mo- 
mentum will be marginal compared to the CENBOL contribution. But as Fq is weakest 
amongst all the components of CENBOL flux, and that F r K is directed towards the axis in 
a larger domain, thus for i-kji c ~ 1 it is possible to observe greater collimation. 

It is to be noted in Figs. (5a-6j), that the various components of radiative moments com- 



puted in thi s paper are quantitatively d ifferent from earlier work s on thin disc (jTajima fc Fukue 



1996 



19981 ). as well as, on slim disc (jWatarai fc Fukue 



J9) 



1999). When comparing radiation 



fields from thin disc and that from TCAF disc model, the first point to be noticed, is the ge- 
ometry of the two disc model is different, secondly the disc motions are different, and thirdly 
the spatial variation of disc intensity (Ico an d Iko) m TCAF model is different from that 
of purely thin disc. Even the radiative contribution from KD of th e TCAF model to variou s 
radiative moments [Figs. f 6a-6i)|, differs from tha t computed by iTajima fc Fukuei (J1998) . 



The reasons are that (a) in lTajima fc Fukuei (jl998f ). the inner radius of the thin disc is 3r g , 
while in this paper the inner radius of KD is x s (~fewxl0r 9 ), so the jet does not 'see' the 
most luminous part of the KD, instead 'sees' the luminous CENBOL, and (b) the shadow 
effect of CENBOL. The shadow effect of CENBOL on jet material is extensively discussed 
in §3.1, we will just point out here that for jet material at z^h s (x — r)/ (x Q — x s ), radiation 
from KD is completely blocked and for z > h s (x — r)/ (x Q — x s ), the jet material only 'sees' 
a fraction of the o uter rim of the KD. Wit h increasing z, the jet 'sees' more and more inner 
part of the KD. In lTaiima fc Fukuei (J1998), such shadow effect was not considered and hence 
the difference in spatial variation of various radiative compone nts between these two pa pers, 
namel y between Figs. (6a-6j) o f this paper and Figs. (2,3,4) of iTaiima fc Fukuei (J1998). 



In 



Watarai &: Fukuei ()1999l ). radiative moments are calculated for slim disc model, for 



three height to disc-radius ratio (H/R), (a) H/R~0.05, (b) H/R~0.4 and (c) H/R~0.56. In 
the present paper, height to radius ratio of the CENBOL is /i s /x s ~0.57 (for x s = 20r g , but 
is always < 0.6 for any high er x s ), so one might be tempted to think that, cases (b) and (c) 



of 



Watarai fc Fukuei (|l999l ) should be similar to Figs. (5a-5j) of this paper. Components of 
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radiative moments are not solely determined by disc geometry, but also on the disc dynamics 



and it s radiative properties, and the CENBOL and the slim disc models of 



Watarai Fnkuf 



19991 ) differ on both these counts. The radial velocity component of the slim disc considered 
was tv~Ci/r 1//2 (=u st in our case) and the azimuthal velocity to be f^,~c 2 /r 1 / 2 (=w^, and c\, 
c 2 depends on advection parameter, viscosity parameter, the ratio of specific heats, etc). The 
radial velocity of the CENBOL has no explicit analytical expression [computed in Appendix 
(A)], and, is not proportional to r -1 / 2 . In the immediate post-shock region, the radial velocity 
(u s t) of matter in CENBOL is much less than the radial velocity of matter in models (b) and 



(c) of 



Watarai fc Fukuei ()1999l ). but close to the black hole it is higher. Even the nature of 



(in CENBOL) and (in slim disc) are different. The higher velocity components close to 
Xi n , enhances the CENBOL intensity close to the black hole, and ultimately enhancing each 
of the radiative moments near the inner edge of the CENBOL compared to that achieved 
in slim disc. 

This implies that the stream line velocity (~10 -4 ) at the base of the jet above the 
CENBOL, will experience greater driving force but as the jet velocity is small there, the jet 
will experience less drag force. The opposite will be true for angular momentum: the drag 
forces will reduce angular momentum more than the jet may gain from the radiation, as Aj n 
is high at the jet base. 

The other major difference in the two models is the difference in radiative property. The 
CENBOL radiation in the rest frame is uniform, while that of the slim disc model falls like 
r -2 . Thus the intensity at the outer edge of CENBOL is comparatively more significant than 
the slim disc case. This makes the radial flux directed inwards in a larger part of the domain 
above the CENBOL, compared to the region above the slim disc. Inspite of the quantitative 
differences between the moments computed, by CENBOL and the sli m disc, the overall 



ualitative similarity between Figs. (5a-5j) of this paper, and Figs. (4) of I Watarai fc Fukue 
1999) is quite evident. 



4 RESULTS 

The results are obtained by integrating Eqs. (6a, 6d, 6g), in the radiation field of the TCAF 
disc given in Eqs (9a-9e). Apart from the disc-parameters supplied to calculate the radiation 
field from the disc, we also supply the injection radius of the jet ?v As the jets are launched 
from the inner surface of the funnel like region of the CENBOL, the injection height (zi) 
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should be just above the CENBOL inner surface. With no loss of generality we take Z{ = 
ri(h s /x s ) + 0.1 in units of r g . We assume that the outflow is made up of purely electron- 
positron pair plasma. We are interested to study the dependence of terminal speed as well as 
the relative spread of jets on disc parameters such as the £ c , rh k , x s , r in . We define r^jz^ as 
the relative spread of the jets, where r^, are the cylindrical radial and axial coordinates 
at which v — n^. If r^j ' z^ < 0.1, we define the jet as well collimated, 0.1 < r^/ 1 z^ < 0.2 as 
fairly collimated, 0.2 < r^/z^ < 0.3 as poorly collimated and so on. 



4.1 Dependence on £ c & m^: 

In Fig. (7a), v s is plotted with log(z), where the CENBOL luminosity is increased from £ c = 
0.2 (long-dashed), £ c = 0.3 (dashed), £ c = 0.4 (solid) in units of L Edd , while the Keplerian 
luminosity corresponds to rhk = 6 (in units of MEdd) an d x s = 20r g . The parameters which 
are kept constant through out the paper are, Vi n = 10 -4 , Aj n = 1.7, Xi n = 2r g and x a = 500r 9 
. The injection radius for the jet is ri n = 3r g in Figs. (7a-b). We see that the streamline 
velocity v s increases with increasing £ c , but the amount of increase in v s for equivalent 
increase in £ c , decreases. 

This is to be expected as it was well documented in Paper-I, that CENBOL radiation 
is a good accelerator. Let us look at Fig. (7b), to see how the corresponding streamlines 
behave. 

In Fig. (7b), the streamlines (z vs r) are plotted for £ c = 0.2 (long-dashed), £ c = 0.3 
(dashed), £ c = 0.4 (solid), while = 6 is kept fixed. Other disc parameters are x s = 20r g 
and ri n = 3r 9 . As the CENBOL luminosity is increased, the streamlines are spreading. 
Though the spreading is decreasing for equivalent increase in £ c , as in the case of streamline 
velocities. Before going into the reason for this let us probe into further details. 

In Fig. (7c), terminal speed is plotted with £ c for rhk = 2 (solid), m k = 5 (dashed) 
rhk — 8 (long-dashed), other disc parameters being r in = 2r g and x s = 20r g . Similar to Fig. 
(7a), we see that increases appreciably with £ c . For lower values of £ c , Voo increases with 
rhk, but decreases with increasing rhk for higher values of £ c . Although the change in 
due to the chan ge in rhk is small compared to the change due to £ c . This phenomenon was 



also reported in 



Paper-ll . As this is rotating jet, we would naturally try to see how the jet A 



behaves at large distances. 

In Fig. (7d), Aoo is plotted with £ c for rhk — 2 (solid), rhk — 5 (dashed) rhk 
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Figure 7. (a) Variation of v s with log(z). Each of the curves represents £ c = 0.2 (long-dashed), £ c = 0.3 (dashed) and £ c = 0.4 
(solid) in units of I/Edd- (b) The streamlines of the jet solution in the z — r plane, corresponding to £ c = 0.2 (long-dashed), 
£ c = 0.3 (dashed) and £ c = 0.4 (solid). For both the figures fa^ = 6, Ti n = 3r g . (c) Variation of terminal speed »oo with 
l c . Each curve represents = 2 (solid), = 5 (dashed), and = 8 (long-dashed), where r;„ = 2r g . (d) Variation 
of terminal specific angular momentum Aoo with £ c . Each curve represents riij, = 2 (solid), rh^ = 5 (dashed), and fn^ = 8 
(long-dashed), where r; n = 2r g . (e) Variation of terminal speed with Each curve represents £ c = 0.1 (long-dashed), 
£ c = 0.25 (dashed) and £ c = 0.4 (solid) in units of Z/Edd- The injection radius of the jet is = 2r g . (f) Variation of terminal 
sp. angular momentum Aoo with mj.. Each curve represents £ c = 0.1 (long-dashed), £ c = 0.25 (dashed) and £ c = 0.4 (solid) in 
units of I/Edd- The injection radius of the jet is ri n = 2r g . For all the figures Vi n = 10 -4 , \i n = 1.7, Xi n = 2r g , x s = 20r g and 
x = 500r 9 . Values of Vi n , Aj n , Xi n and x is kept constant through out the paper. 
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(long-dashed), other disc parameters being r in = 2r g and x s = 20r g . We see that generally 
Aqo increases with £ c . Within the funnel shape region of the CENBOL, the radiation field 
produces strong drag terms along <ft direction which removes A more than A increased by f^, 
but further away from black hole A reduces so much that the drag terms [ocA(r, z)\ become 
ineffective, and the jet gains angular momentum. So as the CENBOL becomes more and 
more luminous the jet gains angular momentum. Furthermore, for higher values of £ c , Aoo 
increases with rhk- This shows that KD radiation do not remove angular momentum any 
better than the CENBOL radiation. 

In Fig. (7e), is plotted with rh k , parame trized by £ c = 0.4 (solid), £ c = 0.25 (dashed) 
and £ c = 0.1 (long-dashed). As was reported in IPaper-1 has a very weak dependence on 
rhfc. In Fig. (7f), A^ is plotted with rh k , parametrized by £ c = 0.4 (solid), £ c = 0.25 (dashed) 
and £ c = 0.1 (long-dashed). It is evident that A^ has a very weak dependence on rhk- 

Let us now see the effect of radiation on spreading the jet. In Figs. (8a-b) we plot r^j 1 z^. 
In Fig. (8a), r^/zoo is plotted with £ c , parametrized for rhk = 2 (solid), rhk = 5 (dashed) 
rhk = 8 (long-dashed) [same parameters as Fig. (7c-d)], while in Fig. (8b), r^/z^ is plotted 
with rhk, parametrized by £ c = 0.4 (solid), £ c = 0.25 (dashed) and 0.1 (long-dashed) [same 
parameters as Fig. (7e-f)]. The most remarkable contrast is that the jets are more collimated 
with increasing KD luminosity, while it tends to spread with increasing CENBOL luminosity. 
Though one must notice that spreading do not increase monotonically with £ c , but tends to 
decrease with equivalent increase in £ c . There are two features which are quite intriguing, 
(i) contrasting nature of CENBOL and KD radiation field in terms of collimation of the jet, 
and (ii) KD radiation seems to have nominal effect on determining Voo and A^, but plays 
relatively a greater role in collimation. 

The reason that radiation from the KD can collimate the jets better than radiations 
from CENBOL, can be understood from Figs. (5-6) and Eq. (6g). We know that rotating 
matter tends to move away from the axis due to centrifugal force. If the angular momentum 
is reduced by the drag forces then the spreading may be arrested. Close to the CENBOL 
surface the radiation field is dominated by the CENBOL radiation and it produces large 
drag forces along (f>. Further out as the jet starts to 'see' KD radiation, A is reduced to the 
extent the drag force becomes marginal and the jet starts to gain some angular momentum. 
From Figs. (5-6), it is evident that the contributions from KD to the total radiation field is 
much less than that of the CENBOL so, angular momentum removed or added to the jet 
by the radiation from KD is marginal [e.g., Fig. (7f)]. From Eq. (6g), it is clear that the 
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Figure 8. (a) Variation of r^/zoo with t c . Each curve represents tiij = 2 (solid), mj, = 5 (dashed), (b) Variation of r^/ Zac 
with mj,. Each curve represents i c = 0.4 (solid), £ c = 0.25 (dashed) and £ c = 0.1 (long-dashed). The curves are drawn for 
x 3 = 20r g . 

spreading of jets depends as much on centrifugal force, as on f r . From the discussion bellow 
Figs. (5-6), we know that, F r c is weakest amongst all the F^s, while F r K <, F K but F r K > F K , 
and is towards the axis in a larger domain. Thus when £ c increases it doesnot help in 
collimation, but as increases the higher negative values of F K , makes f r more and more 
negative, and so it pushes the jet towards the axis of symmetry and helps in collimation, 
inspite of not removing angular momentum to the extent as CENBOL radiation does. 

4.2 Dependence on rf. 

In Fig. 9, we show the effect of injection radius of the jets. In Fig. (9a), streamlines are 
plotted for = 2r g (solid), r\ = 3r g (dashed), r\ = 4r g (long-dashed), r\ = 5r g (dashed- 
dotted), and r\ = 6r g (long dashed-dotted), where £ c = 0.4, m k = 6 and x s = 20r g . In 
Fig. (9b), corresponding A distribution is shown along the streamlines for r\ = 2r g (solid), 
r\ = 3r g (dashed), r\ = 4r g (long-dashed), r\ = 5r g (dashed-dotted), and r\ = 6r g (long 
dashed-dotted), where £ c = 0.4, thk = 6 and x s = 20r g . From both the figures we see that 
the injection height Z{ changes with as has been discussed at the start of this section. It 
is evident that the as the injection radius is varied, the angular momentum of the jets are 



Radiatively driven rotating pair-plasma jets from two component accretion flows 25 




(b) 



10 12 14 



(e) 



0.0004 



0.0003 



0.0002 





,111,111 


111,111,111 


















: 










y 

/ / / 






/// 




- /// 






if 











i ■ 



1 1 r 



I 
\ 

, \ 
\ \ 

\\\ 

\ \ \ 
\ \ 



\\\\ 
\ \\\\ 
\ \ \ \ 



\ \ V 



log(z) 



2 4 6 



10 12 14 




Figure 9. (a) Streamlines for rt = 2r g (solid), r, = 3r g (dashed), r-j = 4r s (long dashed), r, = 5r 9 (dashed-dottcd) & r; = 6r 9 
(long dashed-dotted); £ c = 0.4, = 6. (b) A vs log(z), corresponding to the previous figure, (c) Dk, vs r,, (d) Aoo vs r-j, (e) 
f^oo vs r», and (f) r^/z oo vs t\; for tH^ — 4 (solid), trip. — 7 (dashed) & tti^ — 10 (long dashed). For all the figures x s — 2,0rg. 
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higher and it spreads further. The angular momentum of the jets at different r$ is same, but 
the radiative moments just above the CENBOL surface decreases with r. In other words the 
jets will gain less A from the radiation field, but at the same time they will lose less A due 
to the drag terms. The net effect is that with increasing r, the jets are of higher angular 
momentum and thus these jets are less collimated. 

On the other hand, jets are generated with very low velocity at the base (~Tj), so the drag 
term is negligible, but as all the radiative moments decreases with r above the CENBOL 
surface, the driving that the jets get due to f r & f z is much less as is increased. In Figs. 
(9c-f), various terminal values are plotted with r iy for = 4 (solid), rhj, = 7 (dashed) 
and rrik = 10 (long-dashed), other parameters being i c = 0.4, x s = 20r 9 . In Fig. (9c), we 
see that decreases with r i; but for fixed values of r i; it increas es with l ower value of m fc 



or in other words higher value of £ c /£k- If has been observed in iPaper-J and also in Fig. 



(7c) of this paper, that for £ c £ 0.22, decreases with increasing m k . In Fig. (9d), on the 
other hand shows that Aqo increasing with increasing rj. In Fig. (9e), the rotational velocity 
-u$i$ j u t u l \ at infinity (t^oo) is plotted with r^. Similar to A^, v^oo also increases with 



v 



2 



A interesting feature is seen if one compares Fig. (9f), with the previous two. In Fig. (9f), 
the relative spread r^jz^ is shown to increase with r i: but the interesting feature is, as one 
increases m k) r^jz^ decreases, but both Aoc as well as increases. This vindicates Eq. 
(6g), which says, f r and f z determines the streamlines along with the centrifugal force. 



4.3 Dependence on x s : 

Until now we have investigated the solutions for x s = 20r g . We now concentrate on the 
dependence of jet solutions on x s . 

In Fig. (10a), t>oo is plotted with x s , for £ c = 0.4 (solid), £ c = 0.25 (dashed) and £ c = 0.1 
(long-dashed), where = 3 and = 2r g are kept fixed. With increasing values of x s , the 
CENBOL intensity decreases and the driving force goes down, producing lesser terminal 
speed. We plot Aoo, corresponding to the cases in Fig. (10a), in Fig. (10b). Similar to v^, 
Aqo also decreases with x s . The reason for this is same as in the previous case, i.e., the 
CENBOL intensity decreases with increasing x s . Interestingly, as l c is increased, increases 
appreciably but A^ is increased marginally. Close to the injection radius, the A of the jet 
decreases rapidly [see, Fig. (9b)], and then when at large distances from the disc, A becomes 
too low so that the drag in the azimuthal direction becomes negligible, the jet starts to gain 
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Figure 10. (a) Variation of foo with x s , £ c = 0.4 (solid), £ c = 0.25 (dashed), £ c = 0.1 (long-dashed), for constant rh k = 3. (b) 
Variation of Aoo with x s , £ c = 0.4 (solid), £ c = 0.25 (dashed), £ c = 0.1 (long-dashed), for constant m k = 3. (c) Variation of Uoo 
with x s , rh k = 2 (solid), m k = 5 (dashed), m k = 8 (long-dashed) for £ c = 0.4. (d) Variation of A,^ with x 3 . m k = 2 (solid), 
m k = 5 (dashed), m k = 8 (long-dashed) for £ c = 0.4. For all the figures = 2r g . 



some angular momentum from the radiation. So increasing £ c will produce higher A^, but 
as the gain in A occurs at distances farther away from the disc, the radiative moments falls 
off anyway, so the increase in Aoo is small. 

In Fig. (10c), Voo is plotted with x s , for rhk = 2 (solid), rhk = 5 (dashed) and rhk = 
8 (long-dashed), with consant l c = 0.4. The general conclusion that i>oo decreases with 
increasing x s is still valid, though increasing rhk decreases v^. Another difference we do notice 
is, the curve of i>oo widens with increasing x s , for higher values of £ c in Fig. (10a), while in 
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Fig. (10c), curves converges with increasing x s . The reason is one and the same, i.e., with 
increasing x s , the KD contribution to all the components of radiative moments decreases 
and the CENBOL contribution dominates. As CENBOL radiation is a good accelerator, so 
for higher x s , increasing £ c produces higher relative to the case where is increased. 

In Fig. (lOd), Aqo is plotted with x s , for = 2 (solid), = 5 (dashed) and fn^ = 8 
(long-dashed), with consant £ c = 0.4. Aoo decreases with x s but has almost no dependence 
on KD radiation. This is because the azimuthal component of the radiative flux produced 
by KD is smaller compared to the other components (see §3.2). 

In Fig. (11a), relative spread r^/zoo is plotted with x s , for the same parameters as Fig. 
(lOa-b), i.e., for £ c = 0.4 (solid), £ c = 0.25 (dashed) and £ c = 0.1 (long-dashed), where 
fn k = 3 and = 2r g are kept fixed. We find r^jz^ decreases with increasing x s . As A^ 
decreases with increasing x s [see, Fig. (10b)], so it is not surprising that there would be 
greater confutation. Still for fixed values of x s , the difference in A^ is marginal, while on 
the other hand the difference in r^jz^ is larger! If we turn our attention to Fig. (lib), 
which is plotted for rhk = 2 (solid), rhk = 5 (dashed) and rhk = 8 (long-dashed), with 
consant £ c = 0.4 [same case as Fig. (10c-d)], it is again seen that, r^/z^ generally decreases 
with x s , at the same time, decrement of r^j z^ for fixed values of x s with increasing rhk is 
larger , although is almost indistinguishable. From Figs. (lOa-d), we know why Voo and 
Aqo decreases with increasing x s . But increasing collimation is partly due to decreasing A^, 
otherwise variation of r^jz^ would have just mirrored the variation of A^. Increasing x s 
makes /£ negative in a larger part of the funnel like region. Simultaneously, though radiative 
contributions by KD becomes weaker, none the less it makes f r K < in a still larger part of 
the domain. Hence the combination of decreasing angular momentum as well as f r < in a 
larger region around the axis, collimates the jets to a greater degree. 

Thus we conclude from Figs. (lOa-d, lla-b), if £ c > 0.2, x s > 20r g then jets with terminal 
properties Voq ^ 0.9c and r^/zoo < 0.1 is possible. 

It will be interesting if these results are contrasted with earlier works. It is natural that 
these results would differ from earlier works, since the radiation field of TCAF disc is dif- 



ferent from either thin (|Taiima &: Fukudll998[ ) or slim disc ((Watarai fc FukueNl999l ). For 



both the disc models thin and slim, radiation field generally spreads the jet, and due to 
radiation drag the r adiation field suppresses the motion along z direction. In another model 



( Fukue et al. 



200 1[ ). where the authors considered inner non-luminous disc and outer lumi- 



nous disc, the jets got collimated with increasing disc luminosity and the angular momentum 
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Figure 11. Variation of rx/Zoo with x s . (a) £ c = 0.4 (solid), l c = 0.25 (dashed), l c = 0.1 (long-dashed), for constant fn^ = 3. 
(b) m,fc = 2 (solid), rhk = 5 (dashed), fn^ = 8 (long-dashed) for i c = 0.4. Variation of Aoo with x s . For all the figures r; = 2r g . 



got decreased too. The reason for collimation is that the injection radius of the jet (ro in 
their nomenclature) is less the inner boundary of the disc Xi n (tb in their nomenclature). 
This makes the radial flux f r < in a larger region above (and bellow) the disc surface, and 
it pushes the jet material towards the axis. Along with this, radiat ion drag reduces angular 
mo mentum, so the iet s are c ollimated. In both the earlier papers iTajima fc Fukuei (J1998) 
and Watarai &: Fukuei (1999), the injection radius is greater than the inner radius of the 
disc, which produces f r >0 which pushes the jets outwards, at the same time the maximum 
of is in and around the injection radius, which manages to spin up the jet further, these 
two effects combines to increase the angular momentum of the jets even more. So the crucial 
point for collimation is whether or not f r < 0. 

TCAF disc model has two radiation sources the CENBOL and the outer KD, and 
r, > Xi n . Jets from TCAF model starts with very low streamline velocity (vi n = 10 -4 ), 
and moderadtely high initial angular momentum (A, n = 1.7= angular momentum of the 
CENBOL). As the radiative energy density (e) and the pressure components are quite 
intense close to the inner edge of the disc, the higher value of Aj n ensures very high drag 
force along 0. As a result, the jet A is reduced appreciably, close to the injection radius [e.g., 
Fig. (9b)]. Apart from reduction of the angular momentum of the jets, the radial flux f r is 
negative, just above the CENBOL surface, and this collimates the jets. On the other hand, 
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as Vi n is very small (i.e. drag along streamline is small), resulting in the jets being accelerated 
very fast, shooting upto 30 — 40r g almost vertically, or sometimes initially pushed towards 
the axis [e.g., Figs. (7b, 9a)]. At above 30 — 40r g above the CENBOL surface f r > 0, and the 
drag along is weakened [as, discussed in §4.1], these facts tend to spread the jet, although 
the spreading is arrested by KD luminosity, as well as, by the weakening of the f@. 

One must also remember that, in jets around TCAF model, > x« n . For ri^Xi n , radia- 
tions from the CENBOL region Xi n — is negligible, and hence f r is more negative near r iy 
but for Ti > Xi n , radiations from the CENBOL region x« n — cannot be neglected, making 
f r is less negative near ?v Thus we see for smaller values of collimation and acceleration 
of jets are better [e.g., Figs. (9a-9f)]. 

Although Aqo increases with £ c , but one must notice that Xoo^Xm [e.g., Figs. (7d, 10b, 
10c)], which points to the fact, that indeed A is reduced by drag force. One must also 
notice, that jets are more collimated by increasing outer KD luminosity [e.g., Fig. (8b)], or 
increasing x s [e.g., Fig, (lib)], basically means that jets are basically collimated by f r < 0, 

(2001). 



as was indicated by 



Fnkne et a!. 



4.4 Dependence on £k/£ c '- 

Until now in this section, we have seen that increasing i c produces faster jets, while increasing 
rrik and x s makes the jets more collimated. We have also seen that increasing produces 
less v oo and the collimation is worse. The information of CENBOL radiation is provided by 
its luminosity, the information of KD radiation is supplied by rrik and x s , so to get a better 
understanding, we now study how the relative proportions of CENBOL and KD luminosity 
affect the jet solutions. 

In Fig. (12a), Voo is plotted with £k/£ c , for £ c = 0.5 (solid), £ c = 0.3 (dashed) and 
t c = 0.1 (long-dashed), and x s = 20r g , = 2r g are kept fixed. As in Fig. (7e), we find 
that the increasing from 10% to 100% of £ c , has a marginal effect on v^. In Fig. (12b), 
corresponding r^/z^ is plotted with £k/£ c , «.e. ; for £ c = 0.5 (solid), £ c = 0.3 (dashed) and 
£ c = 0.1 (long-dashed), and x s = 20r g , Ti = 2r g are kept fixed. Though is marginally de- 
pendent on £k/£ c , foo/ Zoo has a stronger dependence on £k/£ c - In other words, KD radiations 
are not strong accelerators but definitely a better collimator. As has been discussed in §4.1, 
the radial flux due to KD is directed towards the axis in a larger domain, and as the ratio 
(■k/^c ~ 0.5 we find highly relativistic and collimated jets. 




Figure 12. Variation of (a) and Too/z^o (b), with @k/£c, for l c = 0.5 (solid), £ c = 0.3 (dashed) and £ c = 0.1 (long- 
dashed); x s = 20r 9 . Variation of Voo (c) and Tao/zoc (d), with i^l^c, for x s = 10r 3 (solid), x 3 = 20r 9 (dashed) and x 3 = 30r 9 
(long-dashed); £ c = 0.4. For all the figures = 2r 9 . 



We have seen in the preceding sub-section that we achieve better collimation with larger 
shock location. In Fig. (12c), we plot with tk/i c , for x s = 10r g (solid), x s = 20r g 
(dashed) and x s = 30r g (long-dashed), where t c = 0.4 and r\ = 2r g are kept fixed. We find 
for x s = 10r g (solid), varies from above 0.94c to just above 0.93c, on the other hand for 
x s = 30r g (long-dashed), t>oo varies slightly and is always a little less than 0.91c. Interestingly 
for smaller values of x s , increasing £k/£ c , or in other words proportionally increasing the KD 
radiation has relatively larger effect (although much smaller than equivalent increase in £ c ). 
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(a) (b) 




Figure 13. Variation of v^, (a) and r^/z^ (b) with x s , for rt = 2r g (solid), r< = 4r 9 (dashed) and r-j = 6r g (long-dashed). 
For both the figures £ c = 0.4 and l k = 0.2. 

For larger x s , the most luminous part of KD (~4r 9 ) is missing so its contribution to the 
components of radiative moments are quite low, and as a result increasing 4 has a negligible 
effect on v^. In Fig. (12d), corresponding r^j z^ is plotted with 4/4, for x s = 10r g (solid), 
x s = 20r g (dashed) and x s = 30r g (long-dashed), where £ c = 0.4 and r\ = 2r g are kept fixed. 
It is also seen that, decrease of r^/zoo with increasing £k/£ c , is more stronger for smaller 
values of x s (solid) than the larger values (long-dashed). 

Thus we see that, if £ c > 0.2, x s > 20r g , 1 > £k/£ c < 0.5, for rj~2r 9 , we get relativistic 
(t>oo > 0.9c) and collimated jets (r^/ ' < 0.1), which means jets will be better collimated 
in intermediate hard states. 

We already know that higher proportions of KD radiation, higher values of x s etc are 
needed for collimation, while £ c accelerates the jets. Still one should not forget about injection 
radius of the jets, because we have seen that increase in can disturb collimation as well as 
produce slower jets. We now study the effect of r iy on jet solution for higher values of £k/£ c 
and increasing x s . 

In Fig. (13a), is plotted with x s for r\ = 2r g (solid), r\ = 4r g (dashed), r\ = 6r g 
(long-dashed), where £ c = 0.4 and 4 = 0.2, i.e., 4/4 = 0.5. The ratio of CENBOL and 
KD luminosity is assumed to be around 0.5, in order to aid collimation. We now see a new 
feature, for lower r\ (solid) decreases monotonically with x s , but for higher r, (long- 
dashed) Voo at first increases up to x s ~15r 9 and then decreases. From Figs. (5a-j), we see 
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that all the radiative moments peaks around ~2r g , just above the inner surface of the 
CENBOL, and also have very strong gradients as one moves along the CENBOL surface 
towards x s . For larger x s the gradients are smoother, but for smaller x s gradients becomes 
stronger. So when r$(= 6r g ) is larger for small x s (long-dashed), then, the radiative moments 
received at is small, because, the presence of strong gradients in all the moments, makes 
all of them peak around 2r g , and sharply decrease at 6r g . As x s increases, for r, = 6r g , 
the gradients are becoming smoother, and = 6r g relatively becomes closer to the axis for 
larger x s . So, as x s increases from 10r g to 15r 9 , the jet material at = 6r g , receives larger 
values of radiative moments, resulting in increase of Voo. Increasing x s further, decreases 
t>oo, because the CENBOL intensity decreases so much that the driving force of radiation 
decreases anyway. 

In Fig. (13b), r^/zoo is plotted with x s for = 2r g (solid), = 4r g (dashed), = 6r g 
(long-dashed), where £ c = 0.4 and £k = 0.2, i.e., £k/£ c = 0.5. The relative spread decreases 
monotonically with x s , as the collimation depends on whether f r < or f r > 0. With 
increasing x s , f r becomes negative (directed towards the axis) in a larger domain, thus 
helping in collimation. Thus we see that for x s > 22r g , r» < 4r 9 , £ c ~0.4 and lhjl c ~ 0.5, we 
have jets with relativistic terminal speed ^ 0.9c and r^/zoo & 0.2, in other words we have 
collimated and relativistic jets. 



5 DISCUSSION AND CONCLUDING REMARKS 



In this paper, we have studied the interaction of radiation with pair dominated jets from 
TCAF disc model. We have ignored the details of the mechanism of production of pair 
dominated jets. High energy photons can produce particle-antiparticle pairs close to the 
inner edge of a disk. It is well known, if the photon energy hv £ 2mc 2 , then an electron- 
positron pair may be created, where h is the Planck's constant, v is photon frequency and 
m is the electron (or positron) mass. If, on the other hand, electron and positron collide 
it will annihilate each other to produce two Gamma-ray photons, a process called pair 
annihilation. Evidently, to produce pair dominated jets pair pro duction has to dominate pair 
annihilation process, as has been theoretically investigated by 



Yamasaki et a,l. 



Mishra fc Melial (jl993h and 



(1999). Observationally, electron-po sitron jets were sugg ested in galactic 



black hole candidate IN 



ova M uscae 1991/GS 112 



( Mirabel fc Rodriguez! Il998l ). in quasar 3C279 (jWardle et al 



-684 (Sunvaev 



et al 11992ft . GRS 1915+105 



1998). Though there is little 
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doubt about the existence of pair dominated jets, radiative acceleration/collimation of such 
jets on the other hand is a different issue altogether. If the pairs are produced to such 
an extent that the jet medium is optically thick then the process discussed in this paper 
fails. For optically thick medium, radiation drag terms are not there, but disc intensity will 
fall off exponentially. We have not considered such details, because considering such details 
involves self-consistent treatment of the inflow-outflow solutions around black holes. The 



present effort confines itself to extend our earlier work (|Paper-j ) to rotating jets. 



The post-shock torus of the TCAF model, produces normal plasma jets and high energy 
photons. In this paper we have considered the electron po sitron jet s were produced within the 
funnel like region of the post-shock torus. In contrast to iPaper-J . we have compute 
different moments of radiation field for axial and off-axial points. Furthermore, in 
we did not consider the Doppler shift of radiation due to the disc motion. The rotational 
motion of matter on the disc surface will generate an azimuthal component of radiative flux. 
The Doppler effect ter m also i nduces non-uniformity in frequency integrated intensity of 



Pacer-!, 



the CENBOL, while in 



Paper-I 



the CENBOL intensity was uniform. This non-uniformity of 
CENBOL intensity resulted in strong gradients in the radiative moments around the inner 
edge of the CENBOL. 

The motion of the Keplerian disc is primarily rotation dominated and its expression is 
known analytically even in general relativity. While the motion of matter in the post-shock 
torus has no simple analytical expression. Since we were not considering inflow-outflow 
solutions self-consistently, we needed to make a proper estimate of the motion of post-shock 
matter. Close to the black hole the infall timescale is much smaller than the viscous time 
scale so the angular momentum of in-falling matter is almost constant near the black hole. 
Thus we assumed wedge flow for the motion of in-falling matter along the inner surface of 
the post shock torus, and solved geodesic equations. The post-shock surface motion may be a 
little over estimated, as pressure gradient terms are ignored. One could have used Paczyhski- 
Wiita potential, but as this potential blows up at the horizon, the geodesic equations are 
solved in general relativistic realm. Another reason for estimating the inflow velocities in 
general relativity is to ensure that the rotational velocities of in-falling matter should tend 
to become zero as one comes closer to the central object. 

The Doppler effect on CENBOL intensity expression [Eq. (8e)], is correct up to first 
order in u/c. We have chopped off the CENBOL at Xi n = 2r g , and the velocity field of the 
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CENBOL is such that at x in < x < x s , 7„~1. At x < Xj n , 7„ sharply increases. Thus taking 
first order correct Doppler term in Eq. (8e) is consistant. 

In this paper, we solved for three equations to find three variables v s , A and r, for which 
three injection values corresponding to the the three variables were supplied. All the other 
parameters like x in , x s , x Q , £ c , rh k , and the disc velocity components u and uk were supplied 
to compute various independent components of radiation stress tensor or, in other words, 
all the moments of radiation field. 

It has been noticed, that the radiation from Keplerian disc has marginal influence in 
determining v s and A, while has a greater role in determining r. As has been explained above, 
collimation depends partly on reducing angular momentum of the jet as well as pushing it 
towards the axis by the radial component of radiative flux. Within the funnel like region 
radiative flux from the Keplerian disc is negative because of the geometry of the TCAF 
model. Furthermore, the azimuthal component of the flux from it is also weak compared to 
the other components, hence the angular momentum gained by Keplerian radiation is small. 
In contrast, the radial flux from the CENBOL is weakest amongst all the other components, 
although, because of the special directiveness of the CENBOL, the radial flux is also towards 
the axis, close to the black hole. However, because of the small size of the CENBOL, at 
^~100r 9 , it approaches a point source and hence the radial flux is positive and spreads the 
jet. But drag terms along r and also the fact that radial flux from CENBOL is weakest 
amongst its all other components, makes the spreading small. The situation dramatically 
changes if the injection radius is increased. As r\ is increased, the radial flux increases and 
is directed more and more away from the axis. On the top of that, the drag terms along 
cf> becomes weaker hence reduction of angular momentum decreases. The net effect is, the 
jet spreads. On the top of that, with increasing r iy the CENBOL intensity goes down so 
the force driving the jets are less, resulting in slower jets. On the other hand, collimation 
is better with increasing x s . Increasing x s , makes the radial flux from CENBOL become 
directed towards the axis in a greater part of the funnel like region, which pushes the jet 
and helps in collimation, although clS X s IS increased the radiation from Keplerian disc has 
lesser influence in determining r. 

In the present paper we have restricted our analysis to non-rotating black holes only. It 
will be interesting to extend this study, to a Kerr black hole, because the higher efficiency 
around a Kerr black hole will produce more intense radiation field, and may produce terminal 
speeds higher than what we have observed in this paper. As the spectrum from TCAF disc 
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around a Kerr black hole is yet to be computed, so the issue of finding the spatial dependence 
of the disc intensity for a Kerr black hole is an open problem as yet. 

We conclude that, if radiative process is the main accelerating and collimating process 
then, 

(i) Electron-positron jets are accelerated to highly relativistic terminal speed, as well as, 
is collimated by the radiations from two component accretion disc model. 

(ii) The space dependent part of the radiative moments from the post-shock region, dom- 
inates the corresponding moments from the Keplerian disc. 

(iii) The CENBOL radiation is the main accelerating agent, but the Keplerian disc radi- 
ation has marginal influence in acceleration. 

(iv) The drag terms in the azimuthal direction is greater than the radiative flux term in 
the same direction, hence the radiation removes angular momentum from the jet, near the 
jet base. 

(v) Collimation is partly brought about by removal of angular momentum, and partly by 
the inward direction of the radial flux. 

(vi) As the radial flux of the Keplerian disc is towards the axis, Keplerian radiation helps 
in collimation. 

(vii) Collimation is better achieved for larger values of shock location and lower values of 
injection radius. 

(viii) As CENBOL radiation is the main accelerating agent and Keplerian disc radiation 
is a good collimator, we conclude that, if this is the main process for acceleration and 
collimation then, highly relativistic and collimated jets should be observed in intermediate 
hard states (£k/^ c ~ 1), and not in extreme hard states (£fc/£ c -Cl). 

(ix) Drawing concrete conclusion, our study shows, if shock in accretion is between 20r g — 
30r 9 , injection radius r\ < 4r g , CENBOL luminosity £ c ^ 0.2LEdd and £fc/£ c ~0.5, then the 
jets will have terminal speed greater than 90% the velocity of light, and the terminal relative 
spread will be less than 20%. 
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APPENDIX A: ESTIMATION OF POST-SHOCK INFLOW VELOCITY 

The CENBOL is assumed conical, i.e., motion of post-shock flow is assumed wedge flow. Let 
us solve the geodesic equation along the radial coordinate r, such that r 2 = x 2 + y 2 where x 
& y defines cylindrical radial and axial coordinated of CENBOL inner surface. The geodesic 
equation for conical inflow around a non-rotating black-hole is; 
3v r — 

*"dJ + T >^ = °" (A1) 
The radial velocity & azimuthal velocity is defined as 



Radiatively driven rotating pair-plasma jets from two component accretion flows 39 
ul = -^~ and u* = £zV>k (A2) 

Denning u = u st / (1 — u 2 ^) 1 ^ 2 , reduces Eq. (A.l) to; 
du \/2 + {u 2 /2-r + l)u 2 



(A3) 



dr r(r — l^aii 

where, 7? = 1/(1 — u 2 ). 

Equation (A. 3) is solved by supplying the value of u s tin = —0.01 at r in = \Jx 2 s + h? s . As 
u st is the 3- velocity of inflow, so u st < 0, i.e., at r — lr g , u st — — 1, — 0. 



